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Abstract. The Polyakov loop is included in the 5 U (2)^ xSU (2)it chiral quark-meson model by considering the 
propagation of the constituent quarks, coupled to the (cr, n) meson multiplet, on the homogeneous background 
of a temporal gauge field, diagonal in color space. The model is solved at finite temperature and quark baryon 
chemical potential both in the chiral limit and for the physical value of the pion mass by using an expansion in 
the number of flavors Nf. Keeping the fermion propagator at its tree-level, a resummation on the pion propagator 
is constructed which resums infinitely many orders in 1/A'/, where 0(1 /Nf) represents the order at which the 
fermions start to contribute in the pion propagator. The influence of the Polyakov loop on the tricritical or the 
critical point in the jiq — T phase diagram is studied for various forms of the Polyakov loop potential. 



1 Introduction 

The low-energy elfective models of the QCD, such as the 
Nambu-Jona-Lasinio (NJL) model and the chiral quark- 
meson model (QM) are based on the global chiral symme- 
try of the QCD. They are very useful to qualitatively un- 
derstand many aspects related to the spontaneous breaking 
of the chiral symmetry and its restoration at finite tempera- 
ture and density. However, the absence of gluonic effective 
degrees of freedom and the lack of color clustering alter 
the reliability of the quantitative thermodynamic predic- 
tions of these models, such as the equation of state or the 
location of the critical end point (CEP) in the Hq-T phase 
diagram. 

Some information on the quark confinement can be in- 
corporated in the effective models through an effective de- 
gree of freedom, the Polyakov loop, which is a good order 
parameter for the deconfinement phase transition in the ab- 
sence of dynamical quarks. The coupling of the Polyakov 
loop to the chiral effective models mimics the effect of 
confinement by statistically suppressing at low tempera- 
ture the contribution of one- and two-quark states relative 
to the three-quark states. This feature makes the Polyakov- 
loop extended effective models more appropriate for the 
description of the low-temperature phase and for quan- 
titative comparison with the thermodynamic observables 
on the lattice [1-3]. Better agreement is expected up to 
T ^ ( 1 .5 -2)7"^. above which the transverse gluonic degrees 
of freedom dominate in thermodynamic quantities, such as 
the pressure, over the longitudinal ones represented by the 
Polyakov loop. 



Some solutions of the PQM model appearing in the lit- 
erature completely disregard quantum effects. The effect of 
including the quantum fluctuation in the PQM model was 
recently studied in [4-6] using functional renormalization 
group methods and also in [7], where it was shown that the 
inclusion of the fluctuations has a significant effect on the 
location of the CEP, which is pushed to higher values of 
Hq. In this contribution we review the results on the Hq-T 
phase diagram obtained in [7] as a result of including dif- 
ferent forms of the effective Polyakov loop potential, as 
compared to those previously obtained in [8] in the chi- 
ral limit of the two flavor QM using the resummation of 
the perturbative series provided by the large-A^/ approxi- 
mation. Starting with the 0-derivable formalism, the ap- 
proximations done to parametrize and solve the model are 
also discussed. 



2 The model in the large-A^/ approximation 

The ingredients for the Polyakov loop extended quark-me- 
son model (PQM) are the (cr, tt) meson multiplet and the 
constituent quark fields ^ coupled to them. These latter 
propagate on the homogeneous background of a temporal 
gauge field. In order to be able to use a large-A^y expan- 
sion we consider Nf constituent quarks and correspond- 
ingly - 1 pions ( = Nf). Performing some rescaling 
with Nf (see [7,8] for details), which assures the finite- 
ness of the tree-level constituent quark mass niq - gv as 
A^^ — > oo, the Lagrangian of the model reads after separat- 
ing the vacuum expectation value vNf of the cr field as : 
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2.1 The grand potential in tlie 0-derivabie 
approximation 

At finite temperature T = 1//3, after analytical continuation 
to imaginary time t — > ;t, Aq — > /A4, the grand partition 
function Z and the grand potential Q{T,hb) of the spatially 
uniform system defined by (1) are defined as 

Z = tr| exp [ - y8(//o(A4) + Hm - tiBQE)^ = e'^" , (2) 

where jUb is the baryon chemical potential. 7/int is the inter- 
acting part of the Hamiltonian constructed from (1). The 
A4-dependent free Hamiltonian H^iAi) for two quark fla- 
vors u and d reads 

Hq(Aa) = 7/0 + d^x \^u\(x)P^luj(x) + id\(x)P^ldj(x^ , 

(3) 

where Ho is the quadratic part of the Hamiltonian at van- 
ishing and i, j denotes color indices. In the so-called Polya- 
kov gauge A4 is diagonal in color space: A4 - diag(^+, <^_, 
-{<p++(p-)). The conserved baryon charge Qb appearing in 

(2) can be expressed in terms of the particle number oper- 

3 

ators as 2b = I 2 A^^^,;, with A^q,i = Nuj + Njj - Naj - N^^ 
(-1 

ande.^. A^^,, = J d^ x^ujui+d.dij. Then, combining Ho(A4) 
and HbQb, the efl'ect of fermions propagating on the con- 
stant background A4, diagonal in the color space is like 
having imaginary chemical potential for color. Following 
Ref. [9], one introduces color-dependent chemical poten- 
tial for fermions 

j"i,2 = - i<l)±, Hi=Hq + ii</>+ + (f>-), (4) 

where fiq = nsl'i is the quark baryon chemical potential. 

3 

Then, introducing the notation = Hq - HiNq one 

!=1 

can write Z as a path integral over the fields, generically 
denoted by ¥ 
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where Qq is the grand potential of the unperturbed system 
with fermions having color-dependent chemical potential 



(6) 



In the 0-derivable approximation of Ref. [10], which 
is also called two-particle irreducible (2PI) approximation, 
the grand potential Q = Q[G^, Go-, G, v, 0, $] is a func- 
tional of the fuU propagators and field expectation values 
of the form : 



/3Q = U(0, 0) + -m^y^ + N^v* - Nhv 
-(N- 1)^ £[lnG;\k) + D;\k)G„{k)] 
- ^J^[]nG-\k) + D-J(k)G^(k)] 
+ tiD.c ^ [in (k) + D-\k)Gik)] + Tzpi , (7) 

where the trace is taken in Dirac and color space. The tree- 
level propagators of the pion, sigma, and constituent quark 
fields are 

iD;lik) = e- mli^, iD-' {k)=H- m„ (8) 

while G„, G,r, and G are the respective full propagators 
in terms of which the set of 2P1 skeleton diagrams de- 
noted by rjpi = r2n[Gn, Go-, G, v, 0, 0] is constructed. To 
0(1/ ^fN) accuracy this is given by 

- Viv|-ftrz,,, J J y5Gik)y5Gik + p)G„ip) 



A 



-^itTD,c f [ 

2^/N JkJp 



G(k)G(k + p)G^(p), (9) 



where the notation n(k) 



= -i I G„' 

Jp 



{p)G„{k + p) was intro- 



duced. The mesonic part of r2p\ contains the 2P1 diagrams 
of the 0{N) model as given in Eq. (2.13) and Figs. 2 and 4 
of [1 1] and also in Eq. (48) and Fig. 2 of [12]. Notice fliat 
the contribution of the fermions goes with fractional pow- 
ers of and intercalates between the leading order (LO) 
and next-to-leading order (NLO) contributions of the pi- 
ons, which go with integer powers of N. 



2.2 Tlie mean-field Polyakov-loop potential 

U(0,0) appearing in (7) represents one of the different 
versions of effective Polyakov-loop potentials defined in 
the literature for pq = 0, where \0\ - \<P\. The simplest 
effective potential is of a Landau type, constructed with 
terms consistent with the Z3 symmetry [13]: 



0''U^,y{0,0) 
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where the temperature-dependent coefficient which makes 
spontaneous symmetry breaking possible is 

b2(T) =ao + ai (^) + «2 (y ) + «3 (y )' - dD 

Tq is the transition temperature of the confinement/decon- 
finement phase transition, in the pure gauge theory Tq = 
270 MeV. The parameters a,, i - 0, . . .,3 and bj, bn deter- 
mined in [14] reproduce the data measured in pure S U(3) 
lattice gauge theory for pressure, and entropy and energy 
densities. When using this potential in either the PNJL or 
the PQM models the minimum of the resulting thermody- 
namic potential is at > 1 for T ^ oo, which is not 
in accordance with the value coming from the definition 

0(x) = <trcL(x)) INc with L(x) = P exp i drA^ir, x) . 

An effective potential for the Polyakov loop inspired 
by a strong-coupling expansion of the lattice QCD action 
was derived in [15]. Using the part coming from the S U (3) 
Haar measure of group integration an effective potential 
was constructed in [16] 

0^ 1/iog(^, 0)=- ^aiT)00 + b{T) hi [l - 600 

+ 4(03 + - 3(00)2] , (12) 

with the temperature-dependent coefficients 

a{T) = ao+ai ( y)+«2 (yf , b{T) = b, (^)' . (13) 

The parameters a,-, i = 0, 1, 2 and b^ reproduce the thermo- 
dynamic quantities in the pure 5 U (3) gauge theory mea- 
sured on the lattice. Since the logarithm in 1J\og{<t>, 0) di- 
verges as 0, — > 1 the use of this effective potential guar- 
antees that 0, — > 1 for T ^ CO. 

A third potential is the one determined in Refs. [15]: 

jSI/pukuC^, 0) = -b [54e-'"'^00 -I- ln(l - 600 

+ 4(03 -I- 03) - 3(00)2)] , (14) 

where the temperature of the deconfinement phase transi- 
tion in pure gauge theory is controlled by a, while b con- 
trols the weight of gluonic effects in the transition. In this 
case, the parameters a = 664 MeV and b - (196.2MeV)3 
are obtained from the requirement of having a first order 
transition at about T = 270 MeV [17,2]. 

It was shown in [17] that there is little difference in 
the pressure calculated from the three effective potentials 
for the Polyakov loop in their validity region up to 7" ^ 
(1.5 - 2)Tc. The presence of dynamical quarks influences 
an effective treatment based on the Polyakov loop which 
in this case is not an exact order parameter. Defining effec- 
tive Polyakov-loop potentials for nonvanishing chemical 
potential when |0| |0| is somewhat ambiguous [2]. Nev- 
ertheless, at yU^ ^0 we use the Z3 symmetric Polyakov- 
loop potentials given above and include the effect of the 
dynamical quarks along the fines of [2], where using renor- 
mafization group arguments the Nf and fiq dependence of 



the T^o parameter of the Polyakov-loop effective potential 
was parametrized as roOu^.A?/) = T^expi-l/iaobijiq))). 
The parameters were chosen to have To(fXq = 0, A'/ = 2) = 
208.64 MeV. When using the Polyakov loop effective po- 
tential given in ( 1 0) and ( 1 2) we will consider in addition to 
To = 270 MeV two more cases, one with a constant value 
To = 208 MeV and the other with the above-mentioned 
ju^-dependent To taken at Nf = 2. 



2.3 The propagator and field equations and the 
approximations made to solve them 

We are interested in the equations for the two-point func- 
tions and the field equations, which are given by the sta- 



tionary conditions 



6G ~ 6G^ ~ 6G~~ ~6v ~ 60~ 6^' ' ^ ' 

In each of these equations the contribution of the fermions 
is kept only at LO in the large-Nf expansion. This contri- 
bution is Oi in the field equations of and 0, 0(1) in 
the equation for the fermion propagator G, and 0(1 / y/N) 
in the remaining equations, that is the field equation of v, 
and the equations of G^r and Go-. 

The second line of (9) does not contribute to any of the 
equations at the order of interest, and the second term on 
the right hand side of (9) contributes only in the equation 
of the sigma propagator 



1 1 Av^ A r 

iG^J(p)^iD-J(p)+ — --jG„ 
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y/N Jk 



G{k)G{k + p). (16) 



In fact. Go- will not appear in any of the remaining five 
equations. Nevertheless, it plays an important role in the 
parametrization of the model, as discussed in Sec. 2.4. 

In the imaginary time formalism of the finite temper- 
ature field theory, which we use for calculation, the four- 
momentum isk = {ioj„, k), where the Matsubara frequen- 
cies are w„ = 2nnT for bosons, while for fermions they 
depend also on the color due to the color-dependent chem- 
ical potential /z, introduced in (4) and are given by w„ = 
(2n -I- \)nT - ifii. The meaning of the integration symbol 
encountered so far is then 



t/3k 



(17) 



The dependence on and of the fermionic trace-log 
term in the grand potential Q and of the quark-pion setting- 
sun in r2Pi results from the fact that, after performing the 
Matsubara sum, the color trace can be expressed in closed 
form in terms of the mean-field (x-independent) Polyakov 
loop and its conjugate 0. The big difference is that while 
in case of the trace-log the result can be expressed in terms 
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of a modified Fermi-Dirac distribution function 



fi(E) = 



and a similar expression for f^{E), but with i-^ and 
Mq ^ ~l-'q^ in case of the setting-sun integral the result does 
not allow for an interpretation in terms of distribution func- 
tions, as shown in the Appendix of [7] between Eqs. (A35) 
and (A37). 

Due to the complexity of the problem, the set of cou- 
pled equations coming from (15) is only solved using some 
approximations described below. 

1 . As a first approximation we disregard the self-con- 
sistent equation for the fermions arising from 6Q/6G = 0, 
that is 

iG-\k) = iD-\k) - ig" r yiG{p)y,GAp - k), (19) 
Jp 

and simply use the tree-level fermion propagator D{k) in 
the remaining five equations. Within this approximation 
the field equation of v and the pion propagator simplify 
considerably. The contribution of the last but one term of 
(9) to the pion propagator breaks up upon working out the 
Dirac structure into the linear combination of a fermionic 
tadpole T{mq) and a bubble integral I{p;mq). Introducing 
the propagator 

Doik) = 7^-^, (20) 



these integrals are defined as 



T(ma) 



I(p; nig) 



Do(q)Do(q + p) 



(21) 



(22) 



In terms of these integrals given between Eqs. (All) and 
(A18) of the Appendix of [7], where the sum over the color 
indices is explicitly done, one obtains 



h = V 



5 h 



,(23) 



r/(/7; w„) . 



(24) 



By making use of the field equation for v in the equation 
for the pion propagator one finds that iG~^{k = 0) = -h/v, 
which means that the Goldstone theorem is fulfilled. Un- 
fortunately, it turns out that this is only accidental, because 
the Ward identity relating the inverse fermion propaga- 
tor and the proper vertex F^^^ = d^F/diffdiffdn" (see e.g. 
Eq. (13.102) of [18]) 



-izTa{75, iG \p)] = V J;^F ^#(0,p, -p), (25) 



is satisfied only with tree-level propagators and vertices. 
The relation above is violated at any order of the perturba- 
tion theory in the large-A^/ approximation, for in view of 
(19) the corrections to the inverse tree-level fermion prop- 
agator are of 0(1), while the corrections to the tree-level 
n - i// - ij/ vertex are suppressed by l/N. The feature is 
probably a shortcoming of our way of implementing the 
large-A^/ scaling in (1), and we could note find a way to 
overcome it. 

2. A further approximation concerns the self-consistent 
pion propagator (24). In [7] several approximations for G^r 
were discussed, here we consider only two of them. 

In the chiral hmit a local approximation is obtained by 
parametrizing the pion propagator as 



p^-M^' 



(26) 



which is then used in all of the equations. is determined 
from = -iG~\{p = 0), which gives the gap-equation 

=m^ + ^{v^ + TFiM))-^^TFimg). (27) 

The subscript F denotes the finite part of the cutoff regu- 
larized integrals defined in Eqs. (21) and (22). In view of 
the field equation (23) in the chiral Umit h = 0, one has 
- 0. We note that due to their self-consistent nature, 
when (27) is solved, a series containing all orders of 1 / "^V 
is in fact resummed. 

In the case of a physical pion mass, in addition to the 
local approximation (26) and (27) for the pion propagator, 
a nonlocal approximation is derived using an expansion to 
(9(1/ Va^) in the expression of the pion propagator (24) ob- 
tained after exploiting the field equation of v (23): 



V^.^^^^^^^.0(i).(28) 



GAP) 



2Nf 



With this form of the pion propagator the field equation for 
V reads 

m^+^ [v" + Tp{M)) + ^^Jp{M, Mg) 

^Tpimg) = -, (29) 

where in this case = h/v and we have introduced the 
integral 

JiM,mg) = -i [ Gliip)p%ip;mg). (30) 
Jp 

Solving this equation for u shows that this approximation 
also resums infinitely many orders in 1 / V^. 

Before proceeding we note that, since our approxima- 
tion is not self-consistent, we do not attempt to do a proper 
renormaUzation by constructing the counterterms which 
absorb the divergences of the integrals. It was shown in [7] 
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that even when a strict expansion in 1 / is performed 
in the pion propagator equation, due to the fact that the 
fermion propagator is unresummed, different subseries of 
the counterterms are needed to cancel the subdivergences 
of different equations. Here we retain the finite parts of the 
integrals obtained using cutoff regularization. We refer the 
interested reader to [7] for details. 



2.4 Parametrization 

We need to determine the mass parameter m^, the cou- 
plings g, A of the Lagrangian (1), the renormalization scale 
M()B, the vacuum expectation value uq = u(T = 0,/ig = 0), 
and the external field h, which vanishes in the chiral limit. 
This is done at T = fiq = using in addition to the pion 
decay constant f„ = 93 MeV, pion mass m„ = 140 MeV, 
and the constituent quark mass taken to be Mg = mN/3 = 
313 MeV, some information coming from the sigma sec- 
tor, such as the mass and the width of the sigma particle 
and the behavior of the spectral function. 

U() is determined from the matrix element of the axial 
vector current between the vacuum state and a one-pion 
state, which due to the rescaling of the vacuum expectation 
value by y/N gives vq - frr/2. The value of the Yukawa 
couphng g = 6.7 is obtained by equating the tree-level 
fermion mass Mq = gvo with Mq. The parameters A and 
Mqb are determined from the sigma propagator, as will be 
detailed below. Having determined them, in the chiral limit 
np- is fixed from the field equation of v, and in the case 
of the physical pion mass np- is determined from the gap 
equation by requiring = m^, and h is obtained from 
the field equation for u, when the local approximation of 
the pion propagator is used, while when the approximation 
(28) for G„ is used, h is fixed by requiring h = ni^vo, and 
np- is determined from the field equation of v. 

In order to fix A and Mqb, one uses in (16) the tree-level 
fermion propagator, the field equation for u (23), and the 
local approximation (26) for the pion propagator, which 
means that in the case of a physical pion mass we neglect 
for simplicity the second term of (29). After working out 
the Dirac structure one obtains the following form for the 
sigma propagator: 



iG-\p) =P--- 



h Ai? 



1 



V 3 \-AIf{p\M)I6 
+ -y^(4ml-p^)lFip;mq). 



(31) 



The integral Ipip', M), obtained using the local approxima- 
tion (26) for the pion propagator with = ml, can be 
found in Eqs. (10) and (11) of [19] with Mq replaced by 
Mqb, while hip^mq) is given in Eqs. (A16)-(A18) of [7]. 

The self-energy has both in the chiral limit M = and 
for M = m„ two cuts along the positive real axis of the 
complex po plane. These are above the thresholds of the 
pion and fermion bubble integrals, which start at - 4M^ 
and p^ = Aniq, respectively. Above these thresholds the re- 
spective pion and fermion bubble integrals have nonvan- 
ishing imaginary parts. We search for poles of the sigma 
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Fig. 1. The A dependence of the real and imaginary parts of the 
complex sigma pole po = - if^/l and of the Landau ghost 
Ml in the chiral limit and for the physical pion mass indicated 
with label /i on the curves. Ml is shown only in this latter 
case, for in the chiral limit there is very little difference. 



propagator analytically continued between the two cuts to 
the second Riemann sheet in the form iG^^ipo = /ce"^, p = 
0) = 0. The pole is parametrized as po = - jTo-/2, with 
the real and imaginary parts corresponding to the mass and 
the half-width of the sigma particle. 

The solution for Mo- and Fg- is shown in Fig. 1 both in 
the chiral hmit {h = m„ = 0) and for the h ¥= case. Sim- 
ilar to the case of the 0(AO model studied in Ref. [19], in 
the chiral limit the value of Mo- is a little smaller and the 
value of Fcr larger than in the h case. Comparing Fig. 1 
with Fig. 2 of Ref. [19] obtained in the OiN) model, that 
is without fermions, the M^-iA) curve moved slightly up- 
ward, but the Fa-iA) curve moved significantly downward, 
which means that in the present case the phenomenologi- 
cally expected value [20] Ma-jFa- ~ 1 cannot be achieved 
for any value of the coupling A. Another difference is that 
for low values of A there are two poles of Go- on the neg- 
ative imaginary axis in contrast to only one such pole in 
the 0(N) model. These poles approach each other as A in- 
creases and after they collide at a given value of A there are 
two complex poles at higher A, one with positive and one 
with negative real part. The imaginary part of the complex 
pole having positive real part is shown in Fig. 1 for the 
renormalization scale Mqb = 885 MeV. As explained in 
the study done in the chiral limit in [8] for lower values of 
the renormalization scale the scale Mi of the lower Lan- 
dau ghost on the imaginary axis comes even closer to Mo- 
and as a result the spectral function of the sigma is heavily 
distorted. In order to avoid this and based on the ratio of 
M^/F^ we have chosen A = 400 and Mqb = 885 MeV. For 
these values M^. = 456 MeV and F^. = 221 MeV in the 
chiral case, while = 474 MeV and Fg. = 152 MeV for 
the case of a physical pion mass. 
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3 The fiq-T Phase diagram 

The thermodynamics is determined by solving the field 
equations, i.e. (23) and the equations giving the depen- 
dence on T and fiq of the two real mean fields and i>, 
which, when the full fermion propagator is replaced by the 
tree-level one have the form: 



dUi0, 0) 



d0 



XSeA d0 d0 



2(f«(m,)- 



Tf(M)) r-r^+ ^ ^ 



d0 



dS/^'^M, nig) dSl^'\M, rriq) \ 



d0 



d0 



d0 
0, (32) 



where Ek - (k^ + w^)2 and M satisfies the gap equation 
(27) or the relation = hi v. The other equation is simi- 
lar to (32), the only difl'erence is that the derivative is taken 
with respect to 0. The integral in (32) is the contribution 
of the fermionic trace-log integral, while the term propor- 
tional with is the contribution of the quark-pion two- 
loop integral in (9) given in Eq. (A35) of [7]. This term 
is disregarded for simplicity when solving the field equa- 
tions for and 0, and only in one case (see the last row 
of Table 2) the complete equation (32) is solved in order 
to estimate the error made by neglecting it in all the other 
cases. 

The tricritical point (TCP) and the critical end point 
(CEP) are identified as the points along the chiral phase 
transition line of the Hq-T phase diagram where a 1 st order 
phase transition turns with decreasing Hq into a 2nd order 
or crossover transition, respectively. In case of a crossover, 
the temperature of the chiral transition is defined as the 
value where the derivative dv/dT has a minimum (inflec- 
tion point of v(T)), while the temperature of the de- 
confinement transition is obtained as the location of the 
maximum in d0/dT. The transition point in the case of 
a 1st order phase transition is estimated by the inflection 
point located between the turning points of the multival- 
ued curve v(fiq) obtained for a given constant temperature. 
Although the precise definition of the 1st order transition 
point is given by that value of the intensive parameter for 
which the two minima of the effective potential are degen- 
erate, we adopt the definition based on the inflection point 
because we compute only the derivatives of the effective 
potential with respect to the fields and propagators. 



3.1 Phase transition in the chiral limit 

In the chiral limit we solve the field equation (23) using 
the local approximation to the pion propagator (26) with 
= and neglect the term proportional with in (32). 
The critical temperature of the chiral transition r^, and the 
pseudocritical temperature of the deconfinement tran- 
sition at vanishing chemical potential, and the location of 
the TCP are summarized in Table I for various forms of the 
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Fig. 2. Phase diagrams obtained in the chiral limit without and 
with the inclusion of the Polyakov loop. The former has lower 
Ttcp and for the latter we used (/iog(<?, <?) with Tq = 208 MeV 
(upper curves) and with Toijig) (middle curves). The deconfine- 
ment transition line is obtained from the inflection point of <P(r). 



Polyakov-loop potential. With the inclusion of the Polya- 
kov loop T^(j^q = 0) and Jtcp increase significantly com- 
pared with the values obtained earlier in [8] without the 
Polyakov loop, but it has little effect on the value of /u^'"^. 
This increase in the value of T^^il-iq - 0) is basically de- 
termined by the value of the parameter To of the Polyakov 
loop potential, while the value of Jtcp shows no signifi- 
cant variation among different cases having the same value 
of To. One can also see, that as explained in [17], the use 
of the polynomial and logarithmic effective potentials for 
the Polyakov loop, that is (10) and (12), drags the value 
of Tj^i/ig = 0) closer to the value of the parameter To than 
the use of f/Fuku(^, 0) given in (14). In this latter case one 
obtains the smallest value for rxcp- 

For To = 270 MeV the deconfinement transition Une 
in the fiq - T phase diagram is above the chiral transition 
line in all three variants of the effective potential for the 
Polyakov loop. When the logarithmic effective potential 
U\og{0, 0) is used either with a constant To = 208 MeV or 
with the yu^,-dependent To proposed in [2] one finds Td < 
T^ at Hq - 0, but at a given value of the chemical potential 
the deconfinement transition Une crosses the chiral transi- 
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139.0 




(60.7,277.0) 


poly 


270 


185.6 


229.0 


(104.5,261.8) 


poly 


208 


168.2 


176.5 


(96.2,263.4) 


log 


270 


191.4 


209.0 


(109.4,261.2) 


log 


208 


167.6 


162.4 


(102.6,261.2) 


log 




167.9 


162.8 


(84.3,266.9) 


Fuku 




176.5 


193.0 


(99.8,262.2) 



Table 1. The (pseudo)critical temperature (Tj) Tj^ of the (decon- 
finement) chiral transition and the at = 0, and the location of 
the TCP in units of MeV obtained in the chiral limit without the 
Polyakov loop (first row) and with the inclusion of the Polyakov 
loop using various effective potentials summarized in Sec. 2.2. 
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tion line and remains above it for higher values of jiq . This 
is shown in Fig. 2, where the deconfinement transition line 
is obtained from the inflection point of ^{T). In contrast to 
the case of constant To, where basically the deconfinement 
transition line is not affected by the increase of Hq, with a 
//^-dependent the deconfinement transition line strongly 
bends, staying close to the chiral line. The two Unes cross 
just above the TCP. 

In the case when To{iJ.g) is used, the lowering of the de- 
confinement transition results in the shrinking of the region 
of the - T plane for which < T < Td, already ob- 
served in Ref. [21]. Since the quantity measuring the quark 
content inside thermally excited particles carrying baryon 
number shows a pronounced change along the chiral phase 
transition line of this region, the region was identified in 
[17] with the so-called quarkyonic phase, a confining state 
made of quarks and is characterized by a high quark num- 
ber density and baryonic (three-quark state) thermal exci- 
tations. 

Comparing our results on the phase diagram to those 
obtained in the chiral limit of the PNJL model one can no- 
tice difi'erences of both qualitative and quantitative nature. 
In the nonlocal PNJL model of Ref. [22] the deconfinement 
phase transition line starts at^q = below the chiral transi- 
tion line both for a polynomial and a logarithmic Polyakov- 
loop eff'ective potential with Tq = 270 MeV, so that the two 
transition Unes cross at finite fiq. In our case this happens 
only for the logarithmic potential with Tq - 208 MeV, as 
can be seen in Fig. 2. In [22,23] the values of Tj^(jdq = 0) 
and Tjcp are much larger than in our case, while the value 
of Hq'^^ is similar to ours. 

3.2 Phase transition with a physical pion mass 

In the approximation (28) for the pion propagator, which 
resum infinitely many orders in 1/ V^, the phase transi- 
tion at r = turns with increasing fiq from a crossover 
type into a first order transition at some value fiq > Mq, 
and in consequence there is a CEP in the //^ - J phase dia- 
gram. The numerical results are summarized in Table 2 for 
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40.7 


(13.5,328.6) 
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(32.9,328.8) 
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22.3 


(30.6,328.8) 


log 


270 


209.7 
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12.0 


(34.5.329.0) 


log 


208 


168.5 


167.1 


*43.0 


(33.0.328.9) 


Fuku 




195.2 


191.3 


21.2 


(31.8,328.8) 


poly 


208 


188.1 


183.1 


21.4 


(32.2,329.0) 



Table 2. The temperatures and Td of the chiral and deconfine- 
ment transitions, the half- width at half maximum of -dv/dT 
at ju^ = (in the case marked with *, due to the asymmetric shape 
of -dv/dT, the full width is given) and the location of the CEP in 
units of MeV obtained using (28) for the pion propagator without 
and with the inclusion of the Polyakov loop. The contribution of 
the quark-pion setting-sun was kept in (32) only for the result of 
the last row. 
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Fig. 3. Phase diagrams obtained for the physical value of the pion 
mass with the inclusion of the Polyakov loop. For the chiral tran- 
sition line which starts at higher T for yu^ = we used Upoiy(<I>, S>) 
with To = 208 MeV and (28) for the pion propagator, for the other 
two phase diagrams we used the local approximation for the pion 
propagator and U\o^{<t>, i>) with To = 208 MeV (middle curves) 
and with ToO"?) (lower curves). The deconfinement transition Une 
is obtained from the inflection point of 0{T). 



various forms of the Polyakov-loop potential reviewed in 
Sec. 2.2. Increasing the temperature jXq decreases and the 
first order chiral restoration becomes a crossover at a much 
lower temperature Tcep than in the chiral case. The inclu- 
sion of the Polyakov loop increases significantly the value 
of TcEP, but, as in the chiral case, it has little effect on the 
value of yu^^^. Neither the choice of the effective potential 
for the Polyakov loop nor the value of To has a significant 
effect on the value of jJ^^^. The result in the last row was 
obtained by keeping in the field equation of the Polyakov 
loop (32) and its conjugate the contribution of the quark- 
pion setting-sun diagram, while in all other cases only the 
contribution of the fermionic trace-log was kept. Compar- 
ing the result in the last row of Table 2 with that of the 
second row obtained using the polynomial Polyakov-loop 
potential, one sees that the error we make by neglecting the 
setting-sun contribution in all other cases is fairly small. 

The values of and 2A jXq = Q are mostly in- 
fluenced by the choice of the Polyakov effective poten- 
tial and the value of Tq : they decrease with the decrease 
of To and by using the logarithmic potential instead of 
the polynomial one. Using the polynomial potential with 
To = 270 MeV the confinement transition line in the Hq — T 
plane is above the chiral transition line. As in the chiral 
case, when a logarithmic potential is used with either a 
fixed value Tq = 208 MeV or with a /.i^-dependent Tq, 
the deconfinement transition line starts at yu^ = below 
the chiral one and the two lines cross at some higher value 
of Hq. This can be seen in Fig. 3. When To(jiq) is used the 
two lines go together until they cross each other just above 
the location of the CEP. This yu^-dependent To gives the 
lowest value of Tcep, similar to the results reported in [24] 
and [6]. Because of the much lower value of the Jcep the 
shrinking of the quarkyonic phase is more pronounced than 
in the chiral case, as the deconfinement transition Unes ap- 
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preaches the jXq axis. This is even more the case here, with 
a physical pion mass, since the deconfinement transition 
is a crossover and as such it happens in a relatively large 
temperature interval. However, the quarkyonic phase does 
not vanish completely as happens in [6], where quantum 
fluctuations are included using functional renormaUzation 
group methods. 

4 Conclusions 

Using the tree-level fermion propagator and some approx- 
imations for the self-consistent pion propagator obtained 
within a large-A'/ expansion, we studied in the S U{2)l x 
SU{2)r chiral quark-meson model, in the chiral Umit and 
for the physical value of the pion mass, the influence of 
the Polyakov loop on the chiral phase transition. When the 
local part of the approximate pion propagator resums in- 
finitely many orders in \INf of fermionic contributions 
it is possible to find a CEP on the chiral phase transition 
line of the fiq - T phase diagram. The inclusion of the 
Polyakov loop potential has a significant effect on Tqep 
and practicaUy no efl'ect on fi^^^ obtained in the original 
chiral quark-meson model, that is which does not contain 
the Polyakov loop. Using the logarithmic form f/iog(?P, ^) 
of the effective potential for the Polyakov loop with pa- 
rameter To = 208 MeV a crossing between the chiral and 
deconfinement transition lines was observed, with the lat- 
ter line starting atfig = slightly below the former one. In 
this case the existence of the quarkyonic phase is possible. 

It was shown in [7] that the result of resumming in 
the pion propagator (9(1/ Va^) fermionic fluctuations ob- 
tained with a strict expansion in 1 / ^^N, while keeping the 
fermion propagator unresummed, the phase transition soft- 
ens to the point that there is no CEP in the Hq-T phase 
diagram within a range < /iq < 500 MeV. For this reason 
it is an interesting question to what extent our results in the 
existence and location of the CEP would be modified by 
the use of the self-consistent propagator for fermions, and 
also by considering the more realistic S U{2>)l x S U{2>)r 
chiral quark-meson model. 



6. T. Herbst, J. Pawlowski, and B.-J. Schaefer, 
arXiv: 1008.0081 [hep-ph]. 

7. G. Marko, Zs. Szep, Phys. Rev. D 82, 065021 (2010). 

8. A. Jakovac, A. Patkos, Zs. Szep, and P. Szepfalusy, 
Phys. Lett. B582, 179 (2004). 

9. C. P. Korthals Altes, R. D. Pisarski, and A. Sinkovics, 
Phys. Rev. D 61, 056007 (2000). 

10. J. M. Luttinger and J.C. Ward, Phys. Rev. 118, 1417 
(1960). 

11. D. Dominici and U. Marini Bettolo Marconi, Phys. 
Lett. B319, 171 (1993). 

12. G. Fejos, A. Paticos, and Zs. Szep, Phys. Rev. D 80, 
025015 (2009). 

13. R. D. Pisarski, Phys. Rev. D 62, 1 1 1501(R) (2000). 

14. C. Ratti, M. A. Thaler, and W. Weise, Phys. Rev. D 73, 
014019 (2006). 

15. K. Fukushima, Phys. Lett. B591, 277 (2004). 

16. C. Ratti, S. RoBner, M. A. Thaler, and W. Weise, Eur. 
Phys. J. C 49, 213 (2007). 

17. K. Fukushima, Phys. Rev. D 77, 114028 (2008). 

18. J. Zinn- Justin, Quantum Field Theory and Critical 
Phenomena, 4th edition, (Oxford University Press, 
Oxford, 2002). 

19. A. Patkos, Zs. Szep, and P. Szepfalusy, Phys. Lett. 
B537, 77 (2002). 

20. I. Caprini, G. Colangelo, and H. Leutwyler, Phys. Rev. 
Lett. 96, 132001 (2006). 

21. H. Abuki, R. Anglani, R. Gatto, G. NardulU, and M. 
Ruggieri, Phys. Rev. D 78, 034034 (2008). 

22. C. Sasaki, B. Friman, and K. Redlich, Phys. Rev. D 
75, 074013 (2007). 

23. P. Costa, C. A. de Sousa, M. C. Ruivo, and H. Hansen, 
Europhys. Lett. 86, 31001 (2009). 

24. M. Ciminale, R. Gatto, N. D. Ippohto, G. NarduUi, and 
M. Ruggieri, Phys. Rev. D 77, 054023 (2008). 



Acknowledgments 

This work is supported by the Hungarian Research Fund 
under Contracts No. T068108 and No. K77534. 



References 

1. W. Weise, Prog. Theor Phys. Suppl. 174, 1 (2008). 

2. B.-J. Schaefer, J. M. Pawlowski, and J. Wambach, 
Phys. Rev. D 76, 074023 (2007). 

3. B.-J. Schaefer, M. Wagner, and J. Wambach, Phys. 
Rev. D 81, 074013(2010). 

4. E. Nakano, B.-J. Schaefer, B. Stokic, B. Friman, K. 
Redlich, Phys. Lett. B682, 401 (2010). 

5. V. Skokov, B. Friman, E. Nakano, K. Redlich, B.-J. 
Schaefer, Phys. Rev. D 82, 034029 (2010). 



